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Measurements of tt production cross-section with DO experiment 

V. Shary for DO coUaboratiorQ 

The recent measurements of the top anti-top quark pair production cross-section in proton an- 
tiproton collisions at = 1.96 TeV in lepton -I- jets and dilepton (including tau lepton) channels 
are presented. These measurements are based on 1 fb~^ of data collected with the DO experiment 
at the Fermilab Tevatron collider. The measured values are compatible with the standard model 
. . . prediction and have the uncertainty close to the uncertainty of the theoretical prediction ~ 10%. 
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^ '. I. INTRODUCTION 

' In the standard model framework the production cross-section of the top quark top antiquark pair (tt) in proton- 
' antiproton collisions could be calculated in the perturbative QCD approach. The most recent calculations are per- 
^ ' formed at the next-to-leading order with the next-to-next-leading order soft correction (NLO-I-NNLO) [l| or at the 
D ] next-to-leading order with the next-to- leading threshold logarithm correction (NLO-I-NLL) 0. The exact value of 
f-^ • the calculated cross-section depends on the top quark mass (jrit) and the choice of parton distribution functions 
^ (PDFs). For example, choosing CTEQ6.6 parametrization 0] leads to a value of 7.39tg j^ pb (for m* = 172 GeV) 
for 111. In [2 for CTEQ6.5 PDFs |j] the cross-section value is calculated to be 7.6i;jjj!^ (scales) toU (PDFs) pb 
(for rrit = 171 GeV). In both cases uncertainties on the calculated value have two main contributions: uncertainties 
^—1 from the choice of PDF and the renormalization and factorization scales, set to be both equal to nit. The scales 
^ ' uncertainty is determined by varying both scales from mt/2 to 2mt. These scales could be varied simultaneously 
, or independently, and different approaches may affect a lot the scale uncertainty Q. Despite some what different 
• approaches, both groups ^ and give compatible values for the tt cross-section, which are also compatible between 
, different sets of PDFs. The uncertainties on these theoretical predictions are estimated to be less than 8%. 
. ■ In this overview we report the tt production cross-section measurements based on the data collected by the DO 
00 ' detector between August 2002 and February 2006 with an integrated luminosity of about 1 fb~^. The description of 
the DO detector can be found elsewhere [5] • In the standard model the top quark decays to a boson and a b quark 
with a probability close to 100%. Here we consider the final state where one W boson decays to a lepton and another 
one to quarks (lepton -I- jets final state, branching ratio 38%) and the final state where both W bosons decay to leptons 
(dilepton final state). For the dilepton final state we consider separately final states with one hadronically decaying 
tau lepton (3.6%) and final states which contain electron and muon only the latter includes r — > e,/i decays (6.5%). 
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II. ELECTRON, MUON DILEPTON FINAL STATES 



In order to maximize the acceptance we consider the following final states: the ones with well identified electron 
or muon (ee, e/i, /i/i), and the final states with one well identified electron or muon and the reconstructed track (e 
+ track , 11+ track). Selection criteria are optimized in each channel separately to yield the best possible precision. 
Possible overlaps between different channels are removed by applying veto cuts. The typical identification criteria for 
an electron require an isolated cluster in the electromagnetic part of the calorimeter with a reconstructed transverse 
momentum pt > 15 GeV. Electron must be matched to a reconstructed track and located in the region < |?7| < 1-1 
or 1.5 < I77I < 2.5, where is a pseudorapidity defined as 77 = —ln{ta.n{9/2)) and 9 is the polar angle with the proton 
beam. Muon is reconstructed as a track in the muon system, matched to a track from the tracking detector, isolated 
from the other objects in the calorimeter and in the tracking detector. Muons are required to be within \ri\ < 2.0 and 
have pt > 15 GcV. Two b quarks from top quark decays reconstructed as jets in the calorimeter using iterative cone 
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Channel 


Observed iV""*' 


BR 


L(pb-^) 


e 


ee 


16 


3.0 


0.01584 


1036 


8.3 % 


efi nj = l 


16 


10.2 


0.03155 


1046 


3.1 % 


e^i nj=2 


32 


6.7 


0.03155 


1046 


12.4 % 




9 


3.6 


0.01571 


1046 


5.1 % 



TABLE I; Observed event yield, expected background, branching ratio, integrated luminosity and efficiency for dilepton ee, 
(1 and 2 jets selections), /i/i channels. 



Channel 


Observed N""'^ 


BR 


L(pb-^) 


£ 


e+track nj = l 


14 


1.58 


0.1066 


1035 


0.25 % 


e+track nj > 2 


8 


1.83 


0.1066 


1035 


1.31 % 


/i+track nj=l 


1 


1.38 


0.1066 


994 


0.18 % 


/i+track nj > 2 


8 


1.36 


0.1066 


994 


1.08 % 



TABLE II: Observed event yield, expected background, branching ratio, integrated luminosity and efficiency for e+track and 
^ + track channels with 1 and 2 jets selections. 

algorithm with a cone radius of 0.5. Jets should be within I77I < 2.5 and typically are required to have pt > 20 GeV, 
where px is corrected for the jet energy scale including the correction for muons from semileptonic h quark decays. In 
all channels except e/x a cut on the missing transverse energy corrected for the jet energy scale and muons pt allows 
to improve signal-to-background ratio due to the presence of neutrinos from W boson decays. The typical cut value 
is 20 - 40 GeV. Additional topological selections exploiting the differences in the signal and background kinematics 
are applied in each channel. Contribution of the main backgrounds, Z boson events decaying to leptons and diboson 
events from WW, WZ, ZZ production are estimated from MC simulation using Alpgen or Pythia generators. Data 
events are used to estimate instrumental backgrounds originating from jets misidentified as electrons or muons from 
the semileptonic b quark decays and events with large missing transverse energy due to the detector resolution effects. 
Table H] summarizes the observed number of events, expected background and efficiency for the ee, e/i, fj,fi channels. 
This allows to measure the tt cross-section with 22% precision: 

a = 6.8tli {stat)toi (syst) ± 0.4 (lumi) pb {mt = 175 GeV) 

In the lepton -|- track channel the signal-to-background ratio is much smaller than in other dilepton channels. In 
order to improve it an additional requirement that at least one jet is identified as a b quark jet (b-tagging) is applied. 
The b-jet identification at DO is based on the neural network, which combines several parameters sensitive to the 
displaced decay vertices of the B hadrons. A typical cut on the neural network output used by the top analysis 
allow to tag b quark jets with an efficiency near 54% and the mistagging rate (probability to tag a light quark jet) 
near 1%. Table [III summarizes the observed number of events, expected background and efficiency for the e -I- track 
and fi + track final states after applying veto selection on ee, e/i, fifj,. The jet multiplicity spectra for tt signal and 
backgrounds events are shown on fig. [TJ Combination of ee, e/i, /i/i, e -I- track and fi + track final states allows to 
increase the precision of the tt cross-section measurement to 19%: 

o- = 6.2 ± 0.9 {stat)to:7 (syst) ± 0.4 (lumi) pb (mt = 175 GeV) 

Systematics uncertainties in dilepton channels come from several sources: jet energy calibration, identification of 
jets, muons, electrons and tracks, trigger efficiency, instrumental background contribution, background normalization, 
b-tagging efficiency (for lepton-|- track only). Many of these sources contribute at the same level and hence there is no 
"main" systematic uncertainty which can be reduce easily. Instead the laborious work on each source of systematics 
is required for further improvement. 

III. DILEPTON FINAL STATES WITH HADRONICALLY DECAYING TAU LEPTON 

Tau lepton is the heaviest lepton and has the strongest coupling to the Higgs boson. Final states with one hadron- 
ically decaying tau lepton contains much more background than other dilepton final states, mainly from multijet and 
W -I- jets events, but it is interesting to study them because they are sensitive to the contribution from the physics 
beyond the standard model. For example, if the charged Higgs boson in the minimal supersymmetric standard model 
scenario has a low enough mass, the t H^b decays will be allowed and will enhance the tau lepton final states. 

The tau lepton is identified as a narrow jet in the calorimeter with a cone size 0.3 matched to one or more tracks. 
The neural network (NN) is used to distinguish tau leptons and jets. The NN input parameters describe energy 
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FIG. 1: Jet multiplicity spectra for dilepton and lepton+track channels. The last bin is inclusive. 



deposition profile in the calorimeter and track-calorimeter correlations. NN is trained on the MC Z tt sample 
and on the background sample from real data. The performance of neural network has been verified with Z — > rr 
data. For more details see [7|. The event selection for er and /ir final states require at least one tau candidate and 
one electron with pT > 15 GeV or muon with pT > 20 GeV, at least two jets with pT > 20 GeV and leading jet 
Pt > 30 GeV, transverse missing energy between 15 and 200 GeV. Further improvement of the signal-to-background 
ratio is reached by the requirement that at least one jet is b-tagged. The main physics background, W^-l-jets and 
Z+jets, is estimated from MC, and multijet background where jets are misidentified as taus is estimated from data. 
The measured cross-section has a precision about 30%, where all tt dilepton and lepton-|-jets final states considered 
as a signal: 

cr = 8.3tl°g istat)t\-j {syst) ± 0.5 (lumi) pb (m* = 175 GeV) 
To measure the cross-section times branching ratio we consider only lepton + tau final state as a signal: 

a{tt)BR{ti IviTVrbb) = 0.19 ± 0.08 {stat) ± 0.07 {syst) ± 0.01 {lumi) pb 
which is compatible with a standard model prediction 0.128 for electron and 0.125 for muon final states. 



IV. LEPTON + JETS FINAL STATES 



In lepton -I- jets final state one W boson decays to a lepton (electron or muon) and another one to jets. In this 
final state two approaches have been used to separate signal from background: kinematic likelihood and b-tagging. 
In both cases, first we define an inclusive sample selected by requiring exactly one isolated electron or muon with 
Pt > 20 GeV and I77I < 1.1 for electron and \ri\ < 2.0 for muons, missing transverse momentum > 20 (for e -I- jets) 
or > 25 GeV (for p. + jets) and at least three jets with px > 20 GeV and \ri\ < 2.5. The leading jet must have 
Pt > 40 GeV, and the lepton pT and the missing transverse energy vector must be separated in azimuth to reject 
background events with mismeasured particles. At this stage the data sample contains only near 20% of tt events. 
The main background originated from W^-l-jets events and multijet events. W + jets events contain real electrons or 
muons coming from W decays. In multijet events one of the jets is misreconstructed as an electron or produces a 
muon because of the semileptonic b quark decays. We determine the multijet background contribution in the selected 
sample by using data samples with the relaxed electron identification or the muon isolation requirements. 

In the kinematic likelihood analysis the further improvement of the signal-to-background ratio is reached by using a 
kinematic determinant built with several variables exploiting the difference in kinematics between tt and backgrounds 
events. For the events with exactly three jets we also use additional requirement that sum of jets transverse momenta 
is less than 120 GeV. Simultaneous fit of e-t-jet and /i-l-jet channels with exactly three jets or with four and more jets 
allows to extract the number of tt and background events. For this we use the discriminant templates from MC for 
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3 jets 


> 4 jets 




3 jets 


> 4 jets 


Ndata 


1760 


626 




Itag >2tags 


Itag >2tags 


Ntt 


245 ± 20 


233 ± 19 


Nu- 


147 ± 12 57 ± 6 


130 ± 10 66 ±7 


Nw+jets. other 


1294 ± 48 


321 ± 30 


Nw+jets 


105 ±5 10 ±1 


16 ±2 2±1 




227 ± 28 


70± 12 


other 


27 ±2 5±1 


8± 1 2± 1 








multijet 


27 ±6 3±2 


6±3 0±2 








Total 


306 ± 14 74 ± 6 


159 ±11 69 ±7 








hline N^ata 


294 76 


179 58 



TABLE III: Sample composition for lepton±jets final state in likelihood analysis (left) and b-tagging analysis (right). Number 
of tt events is calculated using the cross-section measured by the likelihood or by the b-tagging analyses correspondingly. 

the tt signal and W+ jets background and from data for the multijet background. The sample compositions for event 
in the third and forth jet bins is presented in tab. IIIII The measured cross-section value has a 15% precision: 

o- = 6.6 ± 0.8 {stat) ± 0.4 (syst) ± 0.4 (lumi) pb {mt = 175 GeV) 

The b-tag analysis requires that at least one jet is b-tagged. We determine the number of multijet events as above 
and the number of events expected from other background sources from the number of background events in the 
inclusive sample multiplied by their probability to be b-tagged. We obtain the b-tagging probability from the MC 
simulation corrected for differences in the efficiencies observed in the simulation and in data. The composition of the 
b-tagged samples is given in tab. IIIII . The jet multiplicity spectra for tt signal and backgrounds events are shown on 
fig-El The calculated cross-section has a precision of 12%: 

cr = 8.1 ± 0.5 (stat) ± 0.7 (syst) ± 0.5 {lumi) pb {nit = 175 GeV) 




Jet Multiplicity Jet l\1ultiplicity 

FIG. 2: Jet multiplicity spectra for leptons ± jets events with one b-tagged jet (left) and with at least two b-tagged jets 
(right). The histograms show (from top to bottom) the contributions from ti production, from backgrounds with leptons, 
mainly W + jets, and from the multijet background. 

The breakdown of systematics uncertainties is shown in tab. IIVI The b-tag analysis systematics is larger than the 
kinematic likelihood one, mainly because of the systematics on the b-tagging efficiency, but b-tag analysis has better 
statistical uncertainty due to the better signal to background separation. Combining two results together allows to 
improve overall precision, because these two results are only partially correlated (statistical correlation factor is 0.31). 
This combined cross-section value has the relative uncertainty 11%. This is the most precise measurement of the tt 
production cross-section at DO 8]. 

(7 = 7.4 ± 0.5 (stat) ± 0.5 (syst) ± 0.5 {lumi) pb {mt = 175 GeV). 



V. CONCLUSION 

The summary of the DO tt cross-section measurements is shown on fig. [3] All measurements are compatible with the 
standard model prediction. The uncertainty on the measurement in lepton+jet channel is now close to the theoretical 
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Source 


b-tag 


Likelihood 


Combined 


Selection efficiency 


0.26 pb 


0.25 pb 


0.25 pb 


Jet energy calibration 


0.30 pb 


0.11 pb 


0.20 pb 


b tagging 


0.48 pb 




0.24 pb 


MC model 


0.29 pb 


0.11 pb 


0.19 pb 


Multijet bckg normalization 


0.06 pb 


0.10 pb 


0.07 pb 


Likelihood fit 




0.15 pb 


0.08 pb 



TABLE IV: Breakdown of systematic uncertainties in lepton+jets channel. 



uncertainty and will be improved with more statistics collected by DO experiment. At the end of Run II DO plans to 
collect near 8 fb~^. With such integrated luminosity the statistical uncertainty in lepton+jets channel will be near 2% 
and the total uncertainty will be limited by systematics and, in particular, by the luminosity uncertainty. Assuming 
the same relative systematics uncertainty as for the current measurement, one can expect to have near 8% precision 
in the cross-section measurement at the end of Run II. Improving systematics uncertainty by a factor of two will 
decrease the total uncertainty to approximately 7%. The further improvement is not possible without changing the 
method of luminosity determination. For example, the possible approach could be to use pp ^ Z events to normalize 
cross-section measurement. This approach has an advantage because all correlated systematics uncertainties will be 
canceled, but it also requires a careful investigation of the correlation between theoretical uncertainties of the Z and 
ti production cross-sections calculation. 

Another interesting measurement with the a statistics of Run II is the cross-section ratio of dilepton and lepton-|-jet 
channels. In this ratio the luminosity uncertainty and all correlated systematics uncertainties, e.g. jet energy scale, jet 
and lepton ID identifications, will be partially canceled. The precision on this ratio will be limited by the statistical 
uncertainty in the dilepton channel and could reach 5% or even less. 

Cross-section measurement is also useful for the determination of the top quark mass using a theoretical dependence 
of the tt production cross-section on the mass. Such kind of measurements can not compete in precision with the 
direct measurements, but they are much less dependent on the definition of the top mass parameter in MC generators. 
The first results of such kind of measurement look quite promising, see e.g. Q. 
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FIG. 3: Summary of the DO measurements of the tt production cross-section. 
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